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In the Ginsburg Landau (GL ) theory developed by the author to study possible supersolid in 
4 He, there are two parameters: v is the repulsive interaction between the normal component and 
the local superfluid mode, g is a periodic changing chemical potential for the local superfluid mode. 
If g = 0, the normal solid has a particle-hole symmetry. If g is negative ( positive ), the normal solid 
is a vacancy ( interstitial ) type normal solid. Taking the P-H symmetric normal solid as a reference 
state, we show that when g is sufficiently negative, the vacancy induced supersolid is the ground 
state. For comparison, we also discuss interstitial induced supersolid. The microscopic origins of v 
and g are discussed. Experimental implications are given. 



1. Introduction A solid can not flow. It breaks a 
continuous translational symmetry into a discrete lattice 
translational symmetry. While a superfluid can flow even 
through narrowest channels without any resistance. It 
breaks a global U(l) phase rotational symmetry. A su- 
persolid is a state which breaks both kind of symmetries. 
The possibility of a supersolid phase in 4 He was theoret- 
ically speculated in 1970 [l[. Over the last 35 years, a 
number of experiments have been designed to search for 
the supersolid state without success. However, recently, 
by using torsional oscillator measurement, a PSU group 
lead by Chan observed a marked 1 ~ 2% superfluid com- 
ponent of solid 4 He at ~ 0.2K above 25bar, both when 
embedded in Vycor glass and in bulk 4 He [2j. Very re- 
cent specific heat measurements found a excessive specific 
heat peak around the putative supersolid onset critical 
temperature ~ 100 mK [3J. The authors suggested that 
the supersolid state of 4 He maybe responsible for the 
superfluid component. The PSU experiments rekindled 
extensive both theoretical 0, @, @, 0] and experimental 
[1, [t| interests in the still controversial supersolid phase 
of 4 He. There are two kinds of complementary theoret- 
ical approaches. The first is the microscopic numerical 
simulation 4|. The second is the phenomenological ap- 
proach 0, [7|. At this moment, despite all the theoretical 
work, there is no consensus at all on the interpretation of 
PSU's experiments. In [7J, the author constructed a phe- 
nomenological Ginsburg Landau ( GL ) theory to study 
all the possible phases and phase transitions in He and 
map out its global phase diagram from a unified frame- 
work. In this paper, by using the GL theory developed in 
0, we study the possible microscopic mechanism lead- 
ing to the supersolid and analyze carefully the conditions 
for the existence of the supersolid. We find that SS-v is 
more likely than SS-i. However, in order to be complete, 
we study both cases on the same footing. It is also con- 
structive to compare SS-v with SS-i even though the SS-i 
is unlikely to be relevant to the Helium 4 system. We 
call vacancies induced supersolid as SS-v, interstitials in- 
duced supersolid as SS-i. When the two kinds of SS show 
different properties, we treat them differently, when they 
share the same properties, we treat them just in the same 



notation SS. 

2. Ginzburg-Landau theory of a supersolid: 

Let's start by briefly reviewing the GL theory of a su- 
persolid in [7J. The density of a normal solid (NS) is 



defined as n(x) = no + Y2g n G e 
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riQ + 5n{x) where 
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n*Q = n_g and G is any non-zero reciprocal lattice vec- 
tor. In a normal liquid (NL), if the static liquid structure 
factor S(k) has its first maximum peak at k n , then near 
k n , S(k) ~ - +c ^i_ k i yi ■ If the liquid-solid transition 
is weakly first order, it is known that the classical free 
energy to describe the NL-NS transition is (loj : 

f« = H\ r G\ n d\ 2 ~ W Yl U G 1 n G 2 n G 3 S G 1 +G 2 + G 3 ,0 
G Gi,G2,G3 
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r n + c(G 2 — k 2 ) 2 is the tuning parameter 
controlled by the pressure or temperature. Note that be- 
cause the instability happens at finite wavevector, Eqn[T] 
is an expansion in terms of small parameter ng alone, 
it is not a gradient expansion ! Note that the average 
density no does not enter in the free energy. The GL 
parameters w and u n may be determined by fitting the 
theoretical predictions with experimental data. 

Of course, the Superfluid (SF) to Normal Liquid (NL) 
transition at finite temperature is the 3d XY transition 
described by: 



U = K\Wip\ 2 + t\ip\ 2 +u\ip\ 



(2) 



where tp is the complex order parameter and t is the tun- 
ing parameter controlled by the temperature or pressure. 
The GL parameters K, t, u may be determined by fitting 
the theoretical predictions with experimental data. 

The coupling between n(x) and ip(x) consistent with 
all the symmetry can be written down as: 



fint = gSn(x)\ip(x)\ 2 + v(Sn(x)) 2 \ip(x)\ 2 + ■ 



(3) 



where Sn(x) = n(x) — no = Y^'g n G e%G Note that 
the average density no does not enter in the interaction. 



Eqn[3] is an expansion in terms of two small parameters 
5n(x) and ip(x). The ■ ■ • include higher odd and even 
powers of Sn(x) which are sub- leading to the g and v 
term. We keep both g and v term in the Eqn[31 because 
the g term changes the sign, while the v term is invari- 
ant under under the Particle-Hole ( PH ) transformation 
5n(x) —>■ —Sn(x), so the sign of g makes a difference ! 
Due to the two competing orders between the solid and 
the superfluid, we expect v to be always positive and is 
an increasing function of the pressure p. On the other 
hand, we can view g as a periodic chemical potential 
with average zero acting on ip. it is easy to see the cou- 
pling is attractive g v < for vacancies, but repulsive 
g t > for interstitials. If g = 0, the C-NS has the P-H 
symmetry, let's call this kind of PH symmetric C-NS as 
NS-PH ( Fig.l ). In general, g ^ 0, so there is no particle- 
hole symmetry in the C-NS, there are also two kinds of 
C-NS: ( 1) vacancy like C-NS ( named NS-v ) where the 
excitation energy of a vacancy is lower than that of an in- 
terstitial. (2) interstitial like C-NS ( named NS-i ) where 
the excitation energy of an interstitial is lower than that 
of a vacancy. We expect that in contrast to v, g is an 
intrinsic parameter of solid Helium 4 which depends on 
the mass of a 4 He atom and the potential between the 
A He atoms, but not sensitive to the pressure p. 

The GL equations I1I2I3I are invariant under both the 
translational symmetry x — > x + a, n(x) — > n(x + 
a),n(G) -> n(G)e iG ' a ,ip(x) -> ip(x + a) and the global 
U(l) symmetry ip — > ipe lB . In a NL, ng = 0, < ip >= 0. 
In a SF, n e = 0, < ip >^ 0. In a NS, n Q ^ 0, < ip > = 0, 
while in a supersolid ( SS ), ng ^ 0, < ip >^ 0. From 
the NL side, one can approach both the NS and the SF. 
Inside the NL, t > 0, ip has a gap, so can be integrated 
out from Eqn[3l we recover the NS-NL transition tuned 
by r^ in EqnfT] Inside the NL < n(x) >= no, the den- 
sity fluctuations of 5n(x) is massive, so can be integrated 
out from Eqn[3l then we recover the NL to SF transition 
tuned by t in Eqn[2l 

3. The effects of the coupling constants g and v : 
Inside the NS, the translational symmetry is already bro- 
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ken, we can simply set Sn(x) = n(x) — no 
and put it into Eqn[3] to look at the effects of g and 
v. Imagine that at a given pressure p > p c \ ~ 25 bar, 
if tuning g — > 0, but keeping v intact, then the nor- 
mal solid becomes an asymptotically P-H symmetric nor- 
mal solid ( NS-PH ) in Fig.l. It is easy to see that 
v term increases the mass t of ip at k = in the NL 
to t NS - PH = t + vJ2'gHG)\ 2 in the NS-PH. Tak- 
ing the temperature scale t = T — Txy as the refer- 
ence scale, then t^s-PH = T + A(p) where A(p) = 
vT,'g l n (^)| 2 - T XY > is the T = gap for the lo- 
cal superfluid mode ip at k — in the NS-PH. Because 
as the pressure p increases, the repulsive interaction v 
also increases, so it is reasonable to assume that A(p) is 
a monotonic increasing function of p ( Fig. la). Because 
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FIG. 1: (a) The gap A(p) of the local fluctuating superfluid 
mode tjj in the P-H symmetric commensurate normal solid 
( NS-PH ) which exists only at g = 0. We assume it is a 
monotonic increasing function of the pressure p. (b) The zero 
temperature phase diagram of g v verse the pressure p. The 
NS-PH only exists at g v — 0. Any g v 7^ will transfer the 
NS-PH into the NS-v. If the experimental is along the path 
I, then it is a direct 1st order SF to NS-v transition, if it is 
along the path II, it is a 1st order SF to SS-v, then a 2nd 
order SS-v to NS-v transition. 



it is a first order transition across p c i, just like the roton 
gap A r > remains finite just before the first order tran- 
sition, A(p) also remains finite just after the first order 
transition, namely, AfjTjjJ > ( Fig. la ). 

Taking this NS-PH as a reference state, we then grad- 
ually turn on g and see how the ground state evolves. In 
the presence of the periodic potential of n(x) lattice, ip 
forms a Bloch wave, the u interaction of ip in EqnfS] cer- 
tainly favors extended Bloch wave over strongly localized 
Wannier state. In principle, a full energy band calcula- 
tion incorporating the interaction u is necessary to get 
the energy bands of ip. Fortunately, qualitatively impor- 
tant physical picture of GL Eqnsfll [T] [3] can be achieved 
without such a detailed energy band calculation. Ne- 
glecting the self interaction u which is not important in 
the extended Bloch states and taking g ^ as a small 
parameter, a perturbative estimate on the eigenenergy 
e M (0) at the origin k = up to the third order g^ is: 
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where the term linear in g vanishes, fj, = v,i stands for 
vacancies and interstitials respectively and G — 2n/a 
with a ~ 3.17A the lattice constant of the solid 4 He, the 
n(G) = e~ G a / 4a can be taken as a Gaussian where a is 
the width of the Gaussian.. 

In reality, g^ maybe large, so we may need to go beyond 
leading orders. For vacancies, g v < 0, without writing 
out the coefficients explicitly, the expansion is: 



e v (0)=t-g 2 v -\g v \ 3 -\g v \ 4 



(5) 
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so all the coefficients have the same sign. Assuming the 
series converges, we can write e v (0) = t — f v (gv) where 
fv{9v) > /«(0) = is a monotonic increasing function of 
g v and likely has no upper bound. 

For interstitials, gi > 0, without writing out the coef- 
ficients explicitly, the expansion is: 



ei(0) =t-gj +gf - gf + 



(0) 



so the coefficients have oscillating sign. Assuming the 
series converges, we can write £^(0) = t — fi(gi) which 
holds for any gt. Because of the oscillating nature of the 
expansion coefficients, it is hard to judge the nature of 
the function of fi(gi) except that /i(0) = 0. 

The different expansion series of f v (gv) and fi(gi) indi- 
cate that quantum fluctuations may favor vacancies over 
interstitials. However, for simplicity, we assume fi(gt) 
is also a monotonic increasing function of gi, so we can 
discuss vacancies and interstitials induced supersolids on 
the same footing. 

4. Vacancies induced supersolid ( SS-v ): The 
mass of ip v was decreased to t^ v — T + A(p) — f v (g v ) = 
T - Tss-v where T ss _ v {p) = f v (g v ) - A(p) '( Fig.2). 
Because f v {gv) is a monotonic increasing function of \g v \ 
and f v (0) = 0, defining a critical value f v {gvc) = A(p c i), 
then when \g v \ < \g vc \, f v (g v ) < A(p c i), the ip v mode 
remains massive, namely < tp v >= 0. The NS-v remains 
to be the ground state even at T = 0( Fig. lb ). 

If \9v\ > \gvc\, then T S s-v(p c i) = fv(gv) ~ A(p c i) is 
raised above the zero temperature, the SS-v state ex- 
ists in the Fig. lb. Tss-v(Pci) is also proportional to the 
superfluid density measured in the experiments. The re- 
sulting solid is an in-commensurate solid with vacancies 
even at T — whose condensation leads to < ip v >^ 0. 
The SS-v state has a lower energy than the NS-v state at 
sufficiently low temperature. As the pressure increases 
to Pc2 , T ss _ v {p c2 ) = f v (g v ) - A( Pc2 ) = ( Fig.lb ). 
Then f v (g v ) = A(p c2 ), so T SS -v(p) = A{ Pc2 ) - A(p) 
which becomes an effective experimental measure of the 
gap A(p) in the Fig. la. In the following, substituting 
the ansatz t/> = ip% + ipz where < ip%(£) >= ae 1 ® 1 and 
< ^ 2 (£) >= e^J2rn=i^me l ^- s where Q m = G are 
the P shortest reciprocal lattice vectors of the lattice 
into EqnlSl we study the effects of n lattice on ip. From 
EqnJ3J we can see n(x) acts as a periodic potential on 
tp. In order to get the lowest energy ground state, we 
must consider the following 4 conditions: (1) because 
any complex ip ( up to a global phase ) will lead to lo- 
cal supercurrents which is costy, so we can take ip to be 
real, so Q m have to be paired as anti-nodal points. P 
has to be even (2) as shown from the Feymann relation 
Q m , m — 1, • • • , P are simply P shortest reciprocal 
lattice vectors, then translational symmetry of the lattice 
dictates that e(K = 0) = e(K = Q m ), ipi and tp2 have 
to condense at the same time (3) The point group sym- 
metry of the lattice dictates A m = A and is real (4) for 



the vacancies case, the interaction is attractive g v < 0, 
from Eqn[3l we can see that the Superfluid Density wave 
( SDW) p — |?A| 2 simple sits on top of the n lattice as 
much as possible. This is expected, because vacancies are 
hopping near the lattice sites, so the Bose condensation of 
vacancies also happen near the lattice sites. From EqnJJl 
the attractive interaction also favors tp(x = R/2) ~ 
where R/2 stands for any interstitial sites which are in 
the middle of lattice points. It turns out that the the 4 
conditions can fix the relative phase and magnitude of 
and i\) 2 to be 02 = Qi> A = a/P, namely: 



2 P/2 

i/j ss ^ v = V>q(1 + cosQ, 

rn — 1 



x) 



(7) 



where "00 = fl e depends on the temperature and pres- 
sure. Note that in contrast to a uniform superfluid, the 
magnitude of ip is changing in space. This field satisfies 
the Bloch theorem with the crystal momentum k = 
and the Fourier components are ip(K = 0) = a, ijj{K = 
Qm) — a/P. They have the same sign and decay in mag- 
nitude. In principle, higher Fourier components may also 
exist, but they decay very rapidly, so can be neglected 
without affecting the physics qualitatively. P = 6, 8, 12 
corresponds to sc, fee, bec lattices respectively [l2} . 

5. Interstitial induced supersolid ( SS-i ): The 
discussion is quite similar to the vacancy case except that 
(1) if gi < gi c , the C-NS is an interstitial-like C-NS ( 
named as NS-i ) where the interstitial excitation energy 
6j is lower than that of the vacancy e v (2) if gi > g lCl 
the resulting solid is an in-commensurate solid with in- 
terstitials even at T = whose condensation leads to 
< ipi 0, the SS-i state exists in the Fig.2. The argu- 
ments to determine the lattice structure of the SS-i goes 
the same as those in the vacancies case except the condi- 
tion (4), for the SS-v, the interaction is attractive g v < 0, 
so the SDW-v simply sits on top of the n lattice. How- 
ever, for the interstitials case, the interaction is repulsive 
gi > which favors tfj(x — 0) ~ 0, so the Superfluid 
Density wave p — \%j:\ 2 can avoid the n lattice as much 
as possible. It turns out that the the 4 conditions can 
fix the relative phase and magnitude of ipi and to be 
#2 = 6i + A = a/P, namely: 
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i>ss-i = -00(1 - p 2^ COsQ r . 
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where ipo — ae l9 depends on the temperature and pres- 
sure. Note that the crucial sign difference from the va- 
cancies case which make the important difference in the 
X-ray scattering between SS-v and SS-i discussed in [ll[ . 

6. Low energy excitations in the SS: So far, we 
only look at the mean field solutions corresponding to 
vacancies and interstitials. Here we discuss excitations 
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FIG. 2: The phase diagram of T versus g at a given pressure 
Pel < p < Pc2- The finite temperature transitions denoted by 
the dashed line II at a given \g v \ > \g vc \ ( or gi > g ic ) is in 
the 3D XY universality class. However, if \g v \ < \g V c\, then 
the ground state at T — is a vacancy-like C-NS( named as 
NS-v ). Similar thing can be defined for an interstitial-like 
C-NS ( named as NS-i ) when gi < gi C . 



above the mean field solutions. It turns out that the ex- 
citations in both cases are the same, so we discuss both 
cases at same time. In the SS-v and SS-i, the BEC wave- 
functions can be written as 



2 ^ 

Ipss = ^o(l ±p2^ COsQ r , 
m — 1 



x), </>o = |^o|e ie (9) 



where ± sign corresponds to SS-v and SS-i respectively. 
Obviously, there are topological defects in the phase 
winding of 9 which are vortices. At T <C Tss, the vor- 
tices can only appear in tightly bound pairs. However, as 
T — > Tg S , the vortices start to become liberated, this pro- 
cess renders < ipo >= above T > Tss- In addition to 
the superfluid 9 mode in the SS states, there are also lat- 
tice phonon modes u in both n sector and tp sector. How- 
ever, it is easy to see that the coupling Eqn[3]is invariant 
under x^x + u, n{G) -► n(G)e ld - a , ip(G) -> V'(G)e i<5 '", 
so the lattice phonon modes in t/j are locked to those 
in the conventional n lattice. This is expected because 
there is only one kind of translational symmetry break- 
ing, therefore only one kind of lattice phonons. If ignor- 
ing the coupling between the superfluid 9 mode and the 
phonon u mode , then the action to describe the NS to 
SS transition in a static lattice is: 



Uo = K NS \Vijjo\ 2 + t NS \ipo\ 2 + u N s\ipo\ 4 



(10) 



where t NS = T - T ss and T S s-v(p) = fv(j!v) ~ = 
A( Pc2 ) - A(p) or Tss-tip) = fifoi) - A(p) = A( Pc2 ) - 
A(p) ( Fig. 2). Obviously, it is still a 3d XY transition 
with a much narrower critical regime than the NL to SF 
transition jjjj]. It was shown in [a] that the coupling to 
the phonon mode u will not change the universality class 
at finite temperature. As the pressure increased to p c2 . 



Tss-v or Tss-i are suppressed to zero ( Fig. lb), the 
system becomes a NS-v or NS-i where < tpo >= 0. The 
zero temperature transition from the SS-v ( SS-i ) to the 
NS-v ( NS-i ) driven by the pressure near p ~ p c2 will be 



studied in llj. A low energy effective action involving 
the superfluid phonon 9, the lattice phonons u and their 
couplings will also be studied in [111 ]. 

7. Conclusions: By using the GL theory developed 
in 0], we analyze the microscopic mechanism and the 
conditions for the existence of a supersolid. There are two 
crucial parameters in the GL theory: v and g in Eqnj3] 
The v is an increasing function of the pressure p and de- 
termines the gap of the NS-PH, while the g is essentially 
a periodically changing chemical potential for the local 
superfluid mode ip. In contrast to v, g is an intrinsic pa- 
rameter of solid Helium 4 which depends on the mass of a 
4 He atom and the potential between the 4 He atoms, but 
not sensitive to the pressure p. Taking the NS-PH state 
as the reference state, we show that a supersolid can be 
realized by adding small number of vacancies to the NS- 
PH. In fact, the temperature Tss-v(p) — A(p c2 ) — A(p) 
determined in the experiments becomes an effective mea- 
sure of the gap A(p) in Fig. la ). In ill]], we will explore 
many important physical consequences due to this single 
coupling. Recently, in a Path integral quantum Monte- 
Carlo calculation, the authors in [5[ concluded that the 
vacancies will simply move to the boundary and phase 
separate. This corresponds to the g lies in the NS-v 
regime in Fig.l and Fig. 2. However, our results do point 
out that there is a large parameter regime in Fig.l and 
Fig. 2 where the vacancy induced supersolid could be the 
favorable ground state. If the intrinsic parameter g does 
fall into this regime, then there is a NS-v to the SS-v 
transition described by Eqn llOl Indeed, very recent spe- 
cific heat measurements found a excessive specific heat 
peak around the putative supersolid onset critical tem- 
perature ~ 100 mK More independent experiments 
are needed to see if g indeed falls into the regime, so the 
SS-v exists in 4 He. 
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